Regularity properties of optimal transportation 
problems arising in hedonic pricing models 



We study a form of optimal transportation surplus functions which 
arise in hedonic pricing models. We derive a formula for the Ma-Trudinger- 
Wang curvature of these functions, yielding necessary and sufficient con- 
ditions for them to satisfy (A3w). We use this to give explicit new exam- 
ples of surplus functions satisfying (A3w), of the form b[x,y) — H{x + y) 
where _ff is a convex function on R" . We also show that the space of equi- 
librium contracts in the hedonic pricing model has the maximal possible 
dimension, a result of potential economic interest. 

1 Introduction 

Fix Borel probability measures fi and on smooth manifolds X and Y, re- 
spectively, and a smooth surplus function b : X x Y R. Monge's optimal 
transportation problem is to find the map F : X ^ Y , pushing fj, forward to 
which maximizes the transportation surplus: 



Assuming some regularity on fi (say absolute continuity with respect to lo- 
cal coordinates), sufficient conditions on b for the existence and uniqueness of 
the minimizer F were found by Levin |25| . building on results of Gangbo and 
McCann [21], Caffarelli [5] and Gangbo [20]. The regularity, or smoothness, of 
F is currently a very hot area of research. Ma, Trudinger and Wang [30] found 
a condition (A3s) on b under which the optimizer must be smooth as long as 
the marginals fi and v are smooth and bounded above and below. Trudinger 
and Wang then proved that a weaker version of this hypothesis, (A3w), is in 
fact sufficient for the regularity of F [31] [3S] and Loeper showed that, even for 
rougher marginals, (A3w) is both necessary and sufficient for the continuity 
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of F [27]. This framework generalizes and unifies a series of earlier regular- 
ity results obtained by Caffarelli [2][3]|3], Urbas [3^, Delanoe [8] [9] and Wang 
|37j . Since then, many interesing results on regularity have been obtained; see, 
for example, [15][13][ll][lH]IlS][Ii[n][IH][Il][I2[IS]I13 An interesting line of 
current research is to find examples of surplus functions satisfying (A3w) and 
(A3s). 

One main goal in the present paper is to determine when (A3w) holds for 
a particular class of surplus functions. Explicitly, we are interested in functions 
of the form: 

b{x,y) ^ sup h{x,z) + g{y,z) (1) 

Our motivation in studying surplus functions of this form comes partially 
from mathematical economics. A recent paper by Chiappori, McCann and 
Nesheim [7] demonstrated that finding an equilibrium in a certain hedonic pric- 
ing model is equivalent to solving an optimal transportation problem with a 
surplus function of this form (see also Ekeland [TT] [T^] and Carlier and Ekeland 
[B] for another approach to this problem). We briefy review this model now. 
Imagine X parameterizes a population of consumer types who are looking to 
buy some specific good (say houses). The different directions in the manifold 
X may represent various characteristics which differentiate among types (for 
example, age, income, family size, etc) and c?/i(x) represents the relative fre- 
quency of types with characteristics x G X. Suppose now that Y parameterizes 
a space of sellers looking to produce and sell the same good (say companies 
looking to build customized houses). Again, the different directions in Y rep- 
resent characteristics differentiating seller types from each other (for example, 
the size of the company and the location of its headquarters) and di'{y) their 
relative frequencies. Suppose that Z represents the space of available goods 
that can potentially be produced (for example, the collection of all houses that 
can feasibly be built, characterized by their size, location, etc.). The space Z 
is often referred to as the space of contracts. The functions h(x, z) and g{y, z) 
represent the preference of consumer x to buy house z and the preference of 
seller y to build house z, respectively. The result of Chiappori, McCann and 
Nesheim implies that the equilibrium coupling of buyers to sellers is also a solu- 
tion to the optimal transportation problem with marginals fi and v and surplus 
h{x, y) given by equation Despite their relevance in mathemtical economics, 
however, optimal transportation problems with surplus functions of this form 
do not seem to have been studied systematically in the literature. 

The (A3w) condition is that a certain tensor - the Ma-Trudinger-Wang 
curvature, defined in the next section - should be nonnegative on a certain 
set of vectors. Our goal here is to find a formula for the Ma-Trudinger-Wang 
curvature of h in terms of g and h. This will then yield necessary and sufficient 
conditions on g and h in order for 6 to satisfy (A3w). 

In particular, our work here will produce sufficient conditions on g and h 
which will ensure that the equilibrium assignment of sellers to buyers in the 
hedonic pricing problem is continuous. In addition, this result should be of 
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independent interest to mathematicians working in optimal transportation, as 
at present there are few known examples of surplus functions satisfying ( A3w) . 
Our result here yields new examples of such functions. In particular, we show 
that for special choices of the preference functions h and g, we obtain b{x, y) = 
H{x + y), where H is a, convex function on R"0 Our work here yields a new 
formula for the Ma-Trudinger-Wang curvature in terms of the Legendre dual 
H* of H] this formula splits into two terms, one of which is always positive. 

In a recent paper by Figalli, Kim and McCann, a variant of (A3w), called 
(B3w), emerged as a central concept in another type of economic problem [T3]. 
This condition also asserts the positivity of Ma-Trudinger-Wang curvature, but 
on a larger set of vectors than (A3w). Let us mention here that, as we derive a 
formula for the Ma-Trudinger-Wang curvature, our method also yields necessary 
and sufficient conditions on g and h so that b satisfies (B3w). 

In the next section, we recall various structural hypotheses of the function 
b which arise in the regularity theory of optimal transportation, (A3w) and 
(A3s) being the most important. In section 3, we derive conditions on g and 
h that ensure b will satisfy these hypotheses. In the final section, we study the 
structure of the set of contracts that gets signed in equilibrium; explicity, we 
prove that the dimension of this set is the same as the dimension of the support 
of the optimal measure. 

It is a pleasure to thank Robert McCann for suggesting this problem to me 
and for useful discussions during the course of this work. 

2 Assumptions and defintions 

We will assume that the domains X, Y and Z can all be smoothly imbedded 
in larger manifolds, such that their closures X, Y and Z are compact. We will 
also make several assumptions on the preference functions g and h. 

1. The smooth manifolds X, Y, and Z all share the same dimension, which 
we will denote by n. 

2. he C^(X X Z) and g e C^(Y x Z). 

3. For each {x,y) the supremum is attained by a unique z(x,y) G Z. 

4. For all {x,y), the nxn matrix D'^^h{x, z{x,y)) + D'^^g{y, z{x,y)) is non- 
singular. 

Economically, z(x, y) is the contract that maximizes the total utility of agents 
X and y; if, in equilibrium, x and y choose to conduct business with each other, 
z{x,y) is the contract they sign. 

Assumption (1) on the dimensions is not entirely necessary. Most of the 
literature on optimal transportation deals with equal dimensional spaces; one 

^Surplus functions of this form are reminiscent of a class of functions studied by Gangbo 
and McCann | 21| . who were interested in minimizing the transportation cost for cost functions 
of the form l{\x — y\), where Z > is strictly concave. 
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exception is a recent paper by the present author on the regularity of optimal 
transportation when the dimensions differ |33| . For economic applications, how- 
ever, it may be desirable to allow these dimensions to differ. The dimensions 
oi X, Y and Z may represent the number of characteristics used in the model 
to differentiate among buyers, sellers and contracts, respectively, and it is cer- 
tainly possible that these will not coincide. For simplicity, we assume here that 
the dimensions are all equal, but note that it is straghtforward to extend the 
analysis to the case dim{X) > dim(Z) > dim(y), using the extensions of the 
conditions (A0)-(A3s) found in [55] . 

Definition 2.1. We say b is (x , y) -twisted if for all x £ X, the mapping y i— 
Dxh{x,y) is injective on Y. We say b is {y , x) -twisted if for all y £ Y , the 
mapping x i— >■ Dyb{x,y) is injective on X. We say b is hi-twisted if it is both 
(x,y) and {y,x) -twisted. 

If b is (x, ?/)-twisted, the map y Dxb{x,y) is invertible on its range and 
we will denote its inverse by &-expj;(-). 

Definition 2.2. We say b is non-degenerate if for all (x, y) the matrix of mixed, 
second order partials D^yb{x, y) is non-singular. 

We will use analogous terminology for g and h; for example, we will say that 
g is (y, 2:)-twisted ii z i-^ Dyg(y, z) is injective. 

The first three regularity conditions formulated by Ma, Trudinger and Wang 
are: 

(AO) b e C^(x X F). 

(Al) b is bi-twisted. 
(A2) b is non-degenerate. 

For sufficiently regular marginals fx and v, (Al) implies the existence and 
uniqueness of a maximizer F, [25 '^U\^T\^. The condition (A2), in turn, im- 
plies that the graph of F is contained in an n-dimensional Lipschitz submanifold 
of the product X xY [5^. 

Our next definition concerns the structure of the domain Y. 

Definition 2.3. We say Y is b-convex if for all x the set Dxb{x, Y) C T*X is 
convex. 

Ma, Trudinger and Wang showed that the 6-convexity of Y is necessary 
for the continuity of the map F for arbitrary smooth marginals /i and v [30] . 
Assuming this condition, as well as, (A0)-(A2) they showed that under, (A3s), 
which we define below, the optimal map F is smooth. Loeper then showed 
that the weakening (A3w) of this condition is necessary and sufficient for the 
continuity of F [27l. 

To formulate (A3w) and (A3s), we will need the following definition. 
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Definition 2.4. Assume (A0)-(A2) hold and that Y is b-convex. Let x G X 
and y E Y . Choose tangent vectors v 6 T^X and u £ TyY . Set q = Dxb{x, y) G 
T*X and p = {D'^yb{x,y)) ■ u G T*X. For any smooth curve /3(s) in X with 

/3(0) = X and ^(0) = v, we define the Ma-Trudinger-Wang curvature of b at x 
and y, in the directions v and u by: 

MTWl^y){v,u) := ^^b{P{s),b-exp^{tp + q)) 

A local coordinates expression for MTW^^^ y^{v,u) was first introduced by 
Ma, Trudinger and Wang [30]. The formulation above is due to Loeper [37], 
who showed that MTW^^ y-^{v,u) is invariant under smooth changes of coordi- 
nates and, when, —b{x, y) is the quadratic cost on a Riemannian manifold, it is 
equal to the sectional curvature along the diagonal. For general smooth surplus 
functions, Kim and McCann [23^ showed that it is the sectional curvature of 
certain null planes, corresponding to a certain pseudo-Riemannian metric. 

We can now state the final regularity conditions: 

(A3vif) For all {x,y) e XxY,v e T^X and u e TyY such that v'^ ■ Dlyb • u = 
we have MTWl"^ y-){v,u) >0. 

(A3s) For all {x,y) e X xY,v e T^X and u e TyY such that ■ Dlyb • u = 
and ^;, u 7^ we have MTW^^ y^{v, u) > 0. 

In subsequent sections, we will often refer to the curve t b-expx[tp + q) &Y 
as a b-segment in Y . We also note here that the conditions (B3w) and (B3s) 
found in [13] are equivalent to (A3w) and (A3s), respectively, without the 
orthogonality condition v'^ ■ D^yb ■ u — 0. 

3 Regularity properties of b 

The aim of this section is to understand when b satsifies certain regularity 
properties, namely (Al), which ensures the existence and uniqueness of the 
optimal map F, (A2), ensuring the rectifiability of the graph of F, 6-convexity 
of Y and (A3w)/(A3s), governing the regularity of the optimal map. 
First, we verify some simple facts. 

Lemma 3.1. The map z{x, y) is continuously differentiable and b is smooth. 



For all X, y we have the following formulae: 

Dxb{x,y) = Dxh{x,z{x,y)) (2) 

Dyb{x,y) = Dyg{y,z{x,y)) (3) 

D,z{x,y) = ~M-\x,y)Dl^h{x,z) (4) 

Dyz{x,y) = -M^^(x,y)Dlyg{y,z) (5) 

Dlyb{x,y) = -DlMx,z)M-\x,y)Dlyg{y,z) (6) 

where M(x, y) := Dl^h{x, z{x, y)) + Dl^g{y, z{x, y)) 
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Proof. Note that, as z(x, y) maximizes z f->- g(x, z 



) + h{y, z), we have 



D^g{x, z{x, y)) + D^h{y, z{x, y)) = 



(7) 



As D^^g{x, z{x, y))+Dl^h(y, z(x, y)) is non-singular by assumption, the Imphcit 
Function Theorem now imphes that z{x,y) is C^. This, in turn, imphes that 
b{x,y) = g{x, z{x,y)) + h(y,z{x,y)) is at least C^. Now note that b{x,y) — 
g{x, z) + h[y, z) > 0, for all z, with equality when z = z{x, y), which implies ([2]) 
and ([3]). Differentiating ([7]) with respect to x and y, respectively, yields Q and 
(O. Differentiating ^ with respect to y and using ^ yields ©■ □ 

We will first prove results about (AO) — (A2). These proofs are based on 
similar arguments, found in j7j and [53]; we recreate them here for the reader's 
convenience. 

Corollary 3.2. If h and g satisfy (AO), then h satisfies (AO) as well. 

Proof. Using ([7]) and the Implicit Function Theorem we find that z{x,y) is 
C^; equations ([2|) and ^ together with the chain rule now imply the desired 
result. □ 

Proposition 3.3. If h and g are non- degenerate, then b is non- degenerate. If 
h is [x , z) -twisted and g is {z,y) -twisted, then b is {x,y) -twisted. 

Proof. The non-degeneracy part of the proposition follows immediately from 



Assume that h is {x,z) twisted and g is {z,y) twisted. Suppose we have 
Dxb{x,yo) = Dxb{x,yi); we need to showj/o = Vi- By ©, we have D^gix, z{x,yo)) = 
Dxg{x, z{x, yi)) and so by the twistedness oi g we have z{x, yo) = z{x, yi). Now, 
for i = 0, 1 we have Dzg{x, z{x, yi)) + Dzh{yi, z{x, yi)) = by [7l Therefore, as 
z{x,yo) = z{x,yi), 

D^h{yQ,z{x,yo)) = -D^g{x, z{x,yo)) = -D^g{x, z{x,yi)) = D^h{yi, z{x,yi)) 

Again using z(a;,yo) = yi), the equality D2/i(?;o, ^(a:^, yo)) = D^h{yi, z{x,yi)) 
and the twistedness of h imply that j/o = TJi desired. □ 

Of course, an analagous result holds if h is (z, a;) -twisted and g is (y,z)- 
twisted and so we immediately obtain the following. 

Corollary 3.4. If h and g satsify (Al), then so does b. If h and g satisfy 
(A2), then so does b. 

Next, we consider the 6-convexity condition. 

Proposition 3.5. Assume Z is h-convex and for all h-segments Zt at x, —Dzh{x, Zt) 
is in the domain of g — expzS')- Then the domain Y is b-convex. 
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Proof. Let Dxb{x,yi) — pi, for i = 0, 1. For all t G [0, 1], we must show that 
there is some yt £ Y such that Dxb{x,yt) = t{pi — po) + po :— pt- Let Zi = 
z{x,yi)] then Dxh{x,Zi) — Dxb{x,yi) = pi for i = 0,1, by [5] in Lemma I^TTl 
The given conditions imply the existence of a zt such that Dxh{x, zt) = pt and 
a yt such that Dzh{x,zt) + Dzg{yt, zt) = 0. Therefore, zt = z{x,yt), and so 
Dxh{x, yt) = Dxh{x, zt) = pt, as desired. □ 

We now provide a formula for the Ma-Trudinger-Wang curvature in terms 
of g and h. 

Theorem 3.6. Let x G X , z E Z . v (z T^X and Zt be a b-segment at x. Set 
yt = g-expzt{-Dzh{x, zt)). Then, for u = yt, MTWj-^ y^^iv, u) is given by 



To prove this result, we will need a simple lemma about the Ma-Trudinger- 
Wang curvature. This lemma is well known, and the resulting formula for the 
Ma-Trudinger-Wang curvature appears in, for example, [22 , but it does not 
seem to be proven explicitly in the literature, so we provide a proof below. 

Lemma 3.7. Let yt be a b-segment at x: Dxb{x,yt) = tp + q. Then, for u — yt , 



V 



'^■{dI^Hx, zt)-Dlzh{x, zt)-{Dlzh(x, zt)+Dlzg{yt, zt)) ^-Dl^h^x, zt)^ • 

(8) 



MTWl,^yJv,u) 



-TP2 '"'^ ■ Dlj){x,yt) -v 
at'' t=o 



Proof. Letting Xs be a curve such that xq — x and xq — v, we have 




^Jv^ ■ DiXx, yt) ■v) + — ^JD.b{-. yt) ■ £o) 

T72 {v'^ ■Dl^Kx,yt)-v) + —^ {{pt + q)-io) 
at'' t=o at'' t=o 

±^ljy^.DlXx,yt)-v) 



{xs'^ ■ Dl^b{xs,yt) ■ Xs + Dxb{xs,yt) ■ x's) 



□ 



We can now prove Proposition 13.61 
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Proof. In light of the preceding lemma, it suffices to show: 

Dlj){x, yt) = Dl^h{x, zt)-Dl^h{x, Zt)-{Dl^h{x, zt)+Dl^g{yt, zt)) ^■Dl^h{x, zt) 
By Lemma [Xn 

D^h{x,y) = Dxh{x,z{x,y)) 

Differentiating this equation with respect to x, noting that z{x,yt) — zt and 
using the formula in Lemma |3. 1 1 for Dxz{x, y) yields the desired result. □ 

Equation ([8]) for the Ma-Trudinger-Wang curvature naturally splits into two 
terms: 

MTWl^y^){v,u)^A + B 

where 



2 

T 



(dIMx,z,)).v, B^-^\^^vJ -{M, 



Here Mt — Dl^h{x, zt) + D^^g{yt, zt), and Vt = Dl.^h{x, zt) ■ v. The first term, 
A, is exactly the Ma-Trudinger-Wang curvature of h. The second term B, can 
be further refined using the chain rule and the formulae for the derivatives of 



the matrix 



Differentiating vf ■ {Mt)~^ ■ vt twice, applying these formulae and using the 
symmetry of Mt and its derivatives yields: 



B = i^-2vt^ ■ M^r^ ■vt'2vt'^ ■ M^r^ ■vt + 4:Vt'^ ■ M^-^MtMf^ ■vt + 

+ vf ■ Mf^MtMf^ ■ Vt - 2vJ ■ Mt^MtM^^MtM^^ ■ vt) |t=o (9) 

Now, note that the maximality oi z h{xt,z) ~\- g{yt,z) at zt — z[x,yt) 
implies that Mj, the second derivative of this map, is negative semi-definite. 
Assumption 4 at the beginning of section 2 asserts that Mt is non-singular, and 
therefore it is negative-definite. The second and last terms in ^ are therefore 
non- negative, due to the negative definiteness of Mt^ and the symmetry of 
Mt^ and Mt- 

As the following examples shows, for appropriate forms of the functions h 
and g, b{x,y) = H{x -f y) becomes an arbitrary convex function of the sum. 
Understanding when functions of this form have positive cross curvature is in- 
teresting in its own right. In this setting, the present approach — essentially 
doing calculations on H* rather than H — has a distinct advantage; 6-segments 
for H correspond to ordinary line segments for the dual variables. That is, 
instead of evaluating H along a 6-segment, we evaluate H* along a line. 
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Example 3.8. Take X = Y = Z = W and set h{x,z) ^ x ■ z - H*{z), 
where H* is the Legendre transform of some smooth, uniformly convex function 
H :R" ^R, andg{y,z) =yz. Thenb{x,y) = H{x + y). Then MTW^^^y^{v,u) 
is given by: 



^ dz^dz^dz^dz^^^^^^^^^ ^ dz^dz^dz^ dz^dz^dz"^ dz^dz^^^ jPaPfc 

i,j,k,l i,j,k^l,a,r 

(10) 

where p — D^H{x + y) ■ u, w — D^H{x + y) ■ v. 

Proof. The curve yt is given by DH[x + yt) = tp + q, or yt — DH*{tp + 
q) — X, and Zt is given by Zt = D^h{x, Zt) = tp + q. Note that A = 0, as 
D^^h{x, Zt) = 0. Turning to B, note that as D'^^h{x, Zt) = I, we have Vt — v, 
and Mt = —D'^H*{tp+q), so that Mt = -D^H*{tp+q)-p, or, in matrix notation, 

{Mt)^j = -Efc afo"'az'' Pk- Similarly, {Mt)^j ^ ~T,k,i o^SrB^PkPi- 

Now, the first three terms in ^ vanish. The final two terms are exactly the 
desired expression. □ 

The formula (jlOp for the Ma-Trudinger-Wang curvature does not seem much 
simpler than the original formula of Ma, Trudinger and Wang [30], but it has 
one important advantage; the second term is always non-negative, because of 
the symmetry of mixed partials and the positive definitness of D^H. We obtain 
immediately, for example, that on any domain where the fourth order derivatives 
of H* are small compared to its third order derivatives and the second order 
derivatves of i/, h satisfies (A3s) (in fact, it satisfies the stronger condition 
(B3s)). 



4 HausdorfF dimension of the space of contracts 

Here we study the structure of the distribution of contracts that get signed in 
equilibrium in the hedonic pricing problem. As in previous sections, we will 
assume for simplicity that dim(A) = dim(y) — dim(Z) :— n, but remark that 
the results may be extended to the case where dim(A) > dim(Z) > dim(F) in 
a straightforward way. 

Assuming that /i assigns zero mass to every set of Hausdorff dimension less 
than or equal to n—1, the twist condition implies the existence of a unique map 
F : X ^ Y solving Monge's optimal transportation problem; as & is C^, this 
map is differentiable almost everywhere. In economic terms, this means that, 
in equilibrium almost every agent x conducts business with a unique agent 
y := F{x); they sign the contract z{x,F{x)). We define the distribution of 
signed contracts nz to be the push forward of /i by the map x i— )■ z{x, F(x)); 
our goal here is to investigate the structure of this measure. Economically, 
the support of this measure can be interpreted as the set of all contracts that 
are executed in equilibrium, while diiz(,z) represents the relative frequency of 
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contracts that are executed in equilibrium. Our result below implies that this 
set has the maximal possible dimension. 

Proposition 4.1. Assume that h is (x, z) -twisted, g is {z ^ y) -twisted and that 
both g and h are non- degenerate. At any point where the map x i— >■ z{x, F{x)) 
is differentiate, it 's derivative has full rank. 

Proof. Fix a point xq where F is well defined and differentiable. Set yo — F(xq) 
and zq = z(xo, yo). It is well known that there is a function m : X — s- R, twice 
differentiable almost everywhere, such that for all x where u is differentiable we 
have 

Du{x) = D^h{x,F{x)) 
Wherever F is differentiable, u is twice differentiable, and we have 

D^u{x) = Dl,b(x, F{x)) + Dlyb{x, F{x))DF{x) (11) 

In particular, the preceding equality holds at x = xq. Set Pq = F)'^'u-{xq) — 
D'^^b{xo,yo); it is well known that Pq > 0. Rearranging pT|) implies 

DF{xo)^{Dlyb{xo,yo)y'Po- (12) 
Now, the derivative of a; i— )■ z{x, F(x)) at xq is 

Dxz{xo, ya) + Dyz{xo, yo)DF{xo) 

Now, using the formulae in Lemma [3.11 together with ([T^ we see that this 
is equal to 

-M~^{xo,yo) ■ Dl^h{xo, zq) + {Dl^h{xo, zq)) ^ ■ Pq 

= M^'^{xQ,yQ) + {Dl^_h{xQ,zo)) ^ ■ Pq ■ {Dl^hixQ, zq)) ^"j ■ Dl^h{xQ, zq) 

(13) 

Now, as M{xQ,yQ) < 0, we have — Af "^(xq, yo) > 0. Also, Pq > 0, hence 
{Dl^h{xQ,ZQ)) ^PQ{Dl^h{xQ,ZQ)) ^ = {Dl^h{xQ,ZQ)) ^ PQ(^{Dl^h{xQ, zq)) > 0. 

Therefore, the term in brackets in (|13p is positive definite and thus invertible, 
which, as Dl^h{xQ, zq) is invertible implies the invertibility of (jl3p . □ 

This proposition immdediately implies that a large set of contracts gets 
signed in equilibrium: 

Corollary 4.1.1. Assuming that the support of fi has Hausdorff dimension n, 
the distribution [iz of signed contracts has n dimensional support. 
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If the dimensions of X, Y, and Z were not the same, but instead we had 
dim(X) > dim(Z) > dim(F), the proof of the proposition would formally be 
the same, but we would have to be careful to interpret the inverses of D'^^h and 
Dy^g as one sided inverses (the modification of non-degeneracy condition found 
in [33] is that these matrices have full rank, so that their one sided inverses 
exist). Of course, in this case, the dimension of the support of would be 
dim(Z). 
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